Glassy dynamics as a melting process 

On melting dynamics and the glass transition, Part II 
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There are deep analogies between the melting dynamics in systems with a first order phase 
transition and the dynamics Irom equilibrium in super-cooled liquids. For a class of Ising spin 
models undergoing a first order transition - namely p-spin models on the so-called Nishimori line - 
it can shown that the melting dynamics can be exactly mapped to the equilibrium dynamics. In this 
mapping the dynamical — or mode-coupling — glass transition corresponds to the spinodal point, 
while the Kauzmann transition corresponds to the first order phase transition itself. Both in mean 
field and finite dimensional models this mapping provides an exact realization of the random first 
order theory scenario for the glass transition. The corresponding glassy phenomenology can then 
be understood in the framework of a standard first order phase transition. 

PACS numbers: 64.70.Q-,75.10.Nr,05.50.-|-q, 



There is still no agreement on what is the fundamen- 
tal principle behind the experimentally observed abrupt 
change in the relaxation time of super-cooled liquids 
when the temperature is lowered Many scenarios 

and theories have been proposed over the time to describe 
the nature of glasses, and over the last few years the the- 
oretical research has been concentrated on the following 
questions. Is there an underlying critical phenomenon 
behind the glass transition or not? Is the glass transition 
a thermodynamic or a purely dynamic notion? What is 
the correct theory of super-cooled liquids? 

The most remarkable experimental fact in the phe- 
nomenology of super-cooled liquids is the extremely fast 
rise of the relaxation time r, that increases by several 
orders of magnitude as the temperature is decreased by 
only a few percents. The relaxation time is tradition- 
ally fitted (and well approximated) by the Vogel-Fulcher- 
Talman law r oc exp [A/{T -Tvft)] B- The extrapo- 
lated temperature Tvft is found to be close the Kauz- 
mann temperature Tx where the extrapolated entropy of 
the super-cooled liquid becomes smaller than the entropy 
of the crystal Q, a fact that has led to speculations on 
the existence of an ideal — yet impossible to observe in 
finite time — glass transition at Tk w Tvft- Another 
fact pointing in the direction of an ideal glass transi- 
tion is the Adam-Gibbs ^ relation r oc exp [C/A5(r)], 
where AS'(r) is the difference between the entropy of the 
crystal and the liquid, this gives a further link between 
thermodynamic and dynamic behaviour. 

In the last decade, a lot of attention has been devoted 
to growing length scales. In most theories the slower- 
than-exponential relaxation time in fragile super-cooled 
liquids follows from the fact that larger and larger regions 
become correlated as the temperature is lowered, so that 
larger ensembles of particles have to be rearranged col- 
lectively to relax the system into equilibrium. However, 
the standard static correlation function does not show 



any sign of such a growing correlation length, and more 
complex correlation functions thus have to be considered. 
Two length scales that are observed to grow significantly 
when the temperature is lowered have been now identi- 
fied. The first is a purely dynamic length associated to 
spacial heterogeneities in the dynamics [^-Q . The second 
length is an equilibrium one associated with the correla- 
tion of a sub-system with its frozen boundaries fToUT^]. 
The existence of diverging correlation lengths points to- 
wards the glass transition being a critical phenomena. 

Is there really a genuine glass transition? Are the time 
and length scales really diverging when approaching the 
glass transition or are they only growing? A definitive 
answer on these questions is difficult to obtain as both 
simulations and experiments are faced with the extremely 
slow dynamics. According to the random first order the- 
ory (RFOT) of the glass transition [T3 - [l8j |. the above 
mentioned time and length scales have a genuine diver- 
gence at the ideal glass transition temperature, although 
the RFOT theory is not free from criticisms (see for in- 
stance [l^ d^l and references therein) . 

In a companion paper [2l| we have discussed that a 
large part of the glassy phenomenology also appears in 
the melting process of a fully ordered phase above an 
ordinary first order phase transition. The bottom-line of 
the analysis is that when the ordered system is brought 
at higher temperature than the melting point, it is inside 
a metastable state from which it needs to escape, and this 
is associated to diverging time and length scales analog 
to those we just discussed for glasses, except that in this 
case the existence of a genuine transition is doubtless. 
There are, however, also important differences between 
melting above an ordinary first order phase transition 
and the glassy dynamics. 

In this work we will show that for a class of spin mod- 
els, both in mean field and finite dimensional systems, 
these differences are washed away and the melting dy- 
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namics can be shown to be exactly equivalent to the equi- 
librium dynamics above the first order phase transition. 
These models are nothing but variants of the Isin g p - 
spin models that have inspired the RFOT theory |23 - t24| , 
and most of our approach is built on the construction of 
Nishimori and collaborators [1^, [26| , in particular Ozeki 
[27L [28l |. Our results are two-folds: (1) we show that 
the standard mean-field approach to the glass transition 
is exactly mappable to a melting problem of some sort; 
(2) we show that there exists a set of finite dimensional 
models where the melting process is equivalent to the 
glassy dynamics in a glass-forming liquid, and that the 
existence of a first-order transition for the melting prob- 
lem implies the existence of a RFOT-like transition for 
the glass. Our results offer an alternative and potentially 
fruitful way of looking at the glass transition problem, all 
from a theoretical, numerical and experimental point of 
view. In particular, many questions about the glass tran- 
sition may be recasted into the more familiar and simpler 
to describe properties of first order phase transitions. 

The paper is organized as follows: In the section I 
we present a class of disordered Ising spin models and 
concentrate on a special line in the temperature/disorder 
plane — the Nishimori line — where the melting problem 
is equivalent to the equilibrium dynamics. In section II 
we discuss the mean-field version of these models. In 
section III we concentrate on a three-dimensional case. 
We summarize and discuss our results in the last section. 



I. WHEN MELTING IS EQUIVALENT TO 
EQUILIBRIUM DYNAMICS 

Following the ideas of Edwards and Anderson, a large 
part of the progress in the theor y of gl asses have orig- 
inated in studies of spin glasses |29l - [3ll |. In particular, 
the Ising p-spin glass 22-24| provides a mean- field the- 
ory for the structural glass transition [T^ - [l7j , which is an 
exact realization of the early landscape picture of Gold- 
stein Here we shall follow this path, although we 
will not restrict ourself to the mean field theory. 

Consider the p-spin Hamiltonian, with p ~ 3 we have 
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(1) 



where the sum is over some triplets of spins (the pre- 
cise details on how the triplets are chosen depends on 
the geometry of the problem: mean-field lattice, finite- 
dimensional grid, etc.), and the interactions Jijk are 
quenched random variables taken from the distribution 

Pp{J^3k) = p5{J,jk - 1) + (f - p)5{J,jk + 1) , (2) 

where p can vary from p = 1/2 (the spin glass case: 
Jijk = ±1 with equal probability) to p = 1 (the ferro- 
magnetic case with Jijk — 1). 

Consider first the pure ferromagnetic model with p — \. 
The ferromagnetic many-body interaction models un- 
dergo usually a first order transition, both in mean field 



[2ll . |2J| and in finite dimension (see the simulations of 
the plaquette models in (33l-[35j|). At high temperature 
the system is in a paramagnetic/liquid state, while for 
low temperature it is in a ferromagnetic/crystal one, and 
a first order transition at Tp separates the two. 

The dynamical behavior of the ferromagnetic many- 
body interaction models, on the other hand, reproduces 
many aspects of super-cooled liquids and their glass tran- 
sition. In particular, crystallization at Tp seems to be 
easily avoidable when cooling down from a large temper- 
ature, and the super-cooled liquid so obtained has all the 
desired glassy phenomenology. There are, however, sev- 
eral problems (actually common to most glassy systems) 
that prevent analytical and conceptual progress: 

(1) The fact that below the ferromagnetic/melting 
temperature Tp the true equilibrium state is given by 
the ferromagnet/crystal makes always any statement on 
equilibrium super-cooled liquid delicate, to say the least. 
Any discussion involving the description of the putative 
ideal glass transition at temperature Tk < Tp is plagued 
by this problem, since there will always be a temperature 
beyond which the nucleation time towards the crystal will 
be larger than the relaxation time in the super-cooled liq- 
uid phase 0, [IHl ■ It would thus be interesting to find a 
model with Tr- = Tp so that the divergence of the equi- 
libration time in the liquid phase could be well defined. 

(2) Glass formers arc known to be hard to simulate, 
since the time to find an equilibrium configuration grows 
faster than exponentially with the temperature. It would 
be really convenient to have an equilibrium configuration 
to start the simulation with at all temperatures. 

(3) Usually in glassy models, it is really difiicult to 
make any statement which is not coming from numeri- 
cal simulations. It would be really convenient to have 
some analytical results and guarantees of genuine phase 
transitions and divergences. 

Quite surprisingly, there is a conceptually simple, and 
rigorous way, to avoid the above difficulties for a class of 
Ising p-spin models if one works on the so-called Nishi- 
mori line and this is the main topic of this paper. 



A. Physics on the Nishimori line 



As discovered by Nishimori [25|, |20| , there is a special 
line, the Nishimori line (NL), in the temperature-disorder 
phase diagram where many results can be established rig- 
orously, on any lattice and thus in any dimension. Fol- 
lowing [2^, [2I] , we start by recognizing that eq. ([2|) can 
be rewritten as 
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with 



e p = 



1-p 



(3) 



The Nishimori line in the plane (T, p) is defined by Kp = 
(3 (see Fig. [1] for an explicit example). 

Consider a given quantity of interest A, say the free 
energy, the energy or the magnetization, that depends 
on the realization of the disorder {Jijfe}, on the inverse 



3 



2 - 

1.5 - 

1 : . 

0.5 - 

— 
0.5 



Liquid 




Glass 



Solid 



0.6 



0.7 O.i 
P 



0.9 



FIG. 1: (color online) Phase diagram of the mean field 3-spin 
Ising spin glass model on a Bethe lattice with coordination 
c = 5 and a fraction of ferromagnetic couplings p. In the 
spin glass case (p = 1/2) — which is the starting point of the 
random first order theory of the glass transition — there is a 
mode-coupling type transition at Tmct = 0.936 (dashed red) 
while the thermodynamic glass transition is at Tk = 0.849 
(full red). In the purely ferromagnetic case (p — 1) the model 
undergoes a first order ferromagnetic phase transition at Tf- 
The ferromagnetic transition extends to lower p (lower blue 
line), as well as the spinodal line of the ferromagnetic phase 
Tspin (upper blue line) . Using a gauge transformation we show 
that the relaxation dynamics for all T > Tk in the spin glass 
case is strictly equivalent to the melting dynamics starting 
from all spins up, Si = 1, on the Nishimori line, eq. Q, (the 
black line crossing the diagram from upper left to lower right). 



temperature /3, and in general also on the spin configu- 
ration {Si}. As is usual for random systems, we want 
to consider the average of this quantity over the disorder 
parameterized by the variable p. 
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(4) 



Following again [2g|, we consider a new set of Ising (±1) 
variables {cTi} and the following gauge transformation 



Si — SiCi ^ 

•hjk — Jijk(7iOj(7k 



(5) 
(6) 



It is easy to check that this leaves the Hamiltonian in- 
variant, i.e. 'H{{Si},{Jijk}) = 'H[{Sia.,},{J.ijk<Ji(Jj(7k})- 
We now apply the gauge transformation and obtain: 



^ 2A^(coshifp)Af 



A{{Jijka.iajak}, {Sia.i}, (3) 



(7) 



where M is the number of interactions. Note that the 
set of values over which is performed the sum X) j ^^-^^ 
y^g. does not change under the transformation either. 
We now average quantity [A\p over all the 2^ possible 



choices of the gauge transformation {ai} to obtain: 



y — 



2^ {cosh KpY' ^j^^ 2*^ 



J2 e-''^''^^^''^i>A{{J,,k'J^a,<Jk},{S^<J^},f3). (8) 

At this point, one recognizes that the denominator of 
the first term is nothing but the annealed average of the 
partition sum for p — 1/2, and the second term is the 
disorder average of Aexp {—KpT-L) for p = 1/2. In other 
words 
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A{{Jijk<Ji(Tjak}, {Siffi}, (3) 



If we are on the Nishimori line K„ 



[Z{Kp)h_ 

(9) 

Cp = /3, we recognize 
the right hand side to be the annealed thermodynamic 
average of the quantity A. If A itself is a gauge invariant 
thermal average of some quantity then from ^ follows 
that on any lattice the quenched thermodynamic average 
on the Nishimori line is equal to the annealed thermody- 
namic average in the fully disordered case (p = 1/2) (72|. 

Let us denote Z{{Jijk}, Kp) the partition function for 
a given realization of disorder p and inverse temperature 

at inverse temperature Kp as 



p. We denote the thermal average with respect to {ai} 



E{..}-4({c7.})e-^p 
Z({J,,fe},ifp) 



(10) 



We now wish to prove the following identity: 

[A{{J.,k}AS^},mp = 

[{A{{Jijk(7tajCrk}, {Sia,}, /3))kJ 



(11) 



It has a simple interpretation (that will get clearer with 
the examples of the next section): consider an instance 
with disorder p and an equilibrium configuration {ai} at 
inverse temperature Kp. Then the disorder average of 
any quantity A{{Jijk}, {Si}, P) is equal to the disorder 
average of A when Jijk Jijk^icrjak and Si SiCTi. To 
prove it consider another gauge transformation 



Si — SiTi , 
C^i — OiTi , 
Jijk — Jijk^i^j^k ■ 



(12) 
(13) 
(14) 



Note that the Hamiltonian, but also JijkaiO'jCrk and SiUi 
are invariant under (|12m4p . Hence applying (|12m4p and 
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averaging over all choices of t; S {±1} we get 
[{A{{J,jkCri'Jj<Tk},{SiCri}, I3))k,] 



E 

{Jijk} 



^{r^} ^ 



2^2*^ (cosh ifp)*=' 

1 ^ Z{{,hjk},Kp) 



{A{{Jijk(Tia.jak}, {Si(Ji},l3))K^ 



2^(coshXp)*^ 



E 

{Jijk} 



2M 



{A{{Jijkcr^ajak}, {Sia^}, I3))kp 



(15) 



Using (0), this proves relation pT|) . 

In particular, on the Nishimori line, i.e. when Kp — /3 
and pnl = 1/(1 + e"^''): one has the following identity 



(16) 



Thus on the Nishimori line all quantities behave the same 
if all spins and all interactions are multiplied by factors 
CTi corresponding to an equilibrium configuration {ui}. 
Note that the above disorder-averaged identities hold for 
any system, in any dimension, with any number of spins. 
Since these disordered systems are self-averaging in the 
thermodynamic limit (both in the mean field [37l.l38l|. and 
in finite-dimension ^]) identity (IT5|) is also valid in the 
thermodynamic limit on the Nishimori line even for one 
given realization of the disorder. 



B. Identities on the Nishimori line 

Let us give several specific examples of the generic 
identities one can obtain on the Nishimori line. Consider 
the equilibrium energy per spin 



e(/3) 



77 E \ ^J^jkSiSJSk 



-m{s}) 



(17) 



Using © and computing explicitly the annealed averages 
one obtains 



M 

[e(/3)]pNL =-ir7tanh(/3) 



N 



(18) 



This was first shown by Nishimori in [2H. Note that 
this is nothing but the average energy of a configuration 
with all spins up. The fact that the energy is an analytic 
function of temperature implies that if a phase transition 
is present along the Nishimori line then it is purely en- 
tropic, i.e. all non-analyticities in the free energy have to 
stem from non-analyticities in the entropy itself. 

Let us now consider the average magnetization. One 
obtains from (jl6l) that 



Pnl 



qsA, (19) 



where the double thermal average {{aiSi)) is both over 
spins {ffi} and {Si}, hence the r.h.s is the equilibrium 
overlap qea (or Edwards- Anderson parameter [31|). The 
equilibrium magnetization on the Nishimori line is equal 
to the average equilibrium overlap (which is the standard 
order parameter in spin glasses) on the Nishimori line, a 
fact again first shown by Nishimori [2^ . 

Another useful identity is obtained when considering 
the so-called Franz-Parisi (FP) potential [13]. It is a use- 
ful tool to understand the properties of glassy states and 
it is defined as the free energy of a system at temperature 
T that has a fixed overlap q with a reference configura- 
tion {ci} that is one of the equilibrium configurations at 
temperature T' 



/fp(/3, /?',?) 



i:K}e-^'^"--»/(/3,{aJ,g) 



(20) 

Under mapping (|lip . we see that the potential at tem- 
perature T with respect to equilibrium at temperature 
T' is simply the free energy of a system with tempera- 
ture T at fixed magnetization with disorder p = pnl (/?')• 
Indeed let A be the free energy at fixed magnetization 



-/3Ar/(/3,m) 



, (21) 



by the mapping (jlip we have 



»7l)]pNL(/3') 



= [/Fp(/^,/3',m)]p„^(;3,) . 



to) e 



-mis,}) 



Pnl(/3') 

(22) 

In particular, ii /3 — /3' the Franz-Parisi potential is noth- 
ing but the equilibrium free energy at fixed magnetiza- 
tion. This shows that the physics relative to the fully 
magnetized configuration is identical to the equilibrium 
physics: this is the particularity of the Nishimori line and 
a crucial point for this paper. 

One can apply the gauge transformation ((SUB]), and the 
identity (IllL t o dynamical quantities as well, as was first 
realized in [23,[2g|. For instance, Glauber, heat-bath or 
any dynamical evolution that satisfies the balance con- 
dition are gauge invariant (this is simply because the 
Hamiltonian itself is invariant, so that the dynamics is 
not affected by the gauge transformation). Consider for 
the sake of the discussion the heat-bath dynamics initial- 
ized in the fully ferromagnetic configuration, Si{0) — 1, 
for a system on the Nishimori line /3 = Kp. The probabil- 
ity that the magnetization evolving under that dynamics 
has a certain value reads 



[P(TO(i)|TO(0) = l)]p,, = 

P{Y.Sm<^,|N\Y,S^{Q)a, 



N] 



Pnl 



Pnl 



.(23) 
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Hence, for the p-spin model on the Nishimori hne, given 
a Hamihonian based dynamics, the decorrelation from 
an equihbrium configuration follows the same functional 
form as the decorrelation from the fully ordered configu- 
ration, a result first shown by [13, HI] . Eq. implies 
that melting dynamics is equivalent to equilibrium dy- 
namics: on the Nishimori line the magnetization decay 
starting from the fully ordered configuration is equal to 
the dynamical correlation function 



(24) 



if Si(ti) = (Ti is an equilibrium configuration. One can 
push this idea one step further and obtain interesting 
identities for the evolution of the dynamic ferromagnetic 
susceptibility xf- Let us denote {-j^jj^jj the average over 
many realizations of the dynamics, then xf reads 

(25) 

Using the same transformation, we see that Xf(^) on the 
Nishimori line equals the equilibrium 4-points suscepti- 
bility X4(0 defined as 

= f E [{^'(^1 + t)S,{t,)S,{h + i)^,(ti)},,„ 

- {S.{h + m{ti)}dy,ASAii + t)Si{h)}dJ (26) 



with Si{ti) = ai being again an equilibrium configura- 
tion. The melting process on the Nishimori line thus 
satisfies ITa 



xTit) = P WW'} 



i dyn 



{m.{t)} 



dyn 



PNL 



xr(t) = /3 {cmayn-{cmi,. 



The bottom line here is that the physics relative to the 
fully magnetized state is at the same time the equilib- 
rium physics so that the melting dynamics from the fully 
ordered state is equivalent to the equilibrium dynamics. 



C. First order transitions on the Nishimori line 

In a companion article (2l| we have discussed analo- 
gies and differences between the equilibrium dynamics in 
super-cooled liquids and the melting dynamics in a sys- 
tem with a general first order phase transition. 

To remind the main findings, we consider the melting 
dynamics in a spin system with a first order phase tran- 
sition. We initialize the system in the fully ordered state 
(that is the completely magnetized one) and suddenly 



change the temperature to put the system in the param- 
agnetic phase. The system will melt into the less ordered 
phase. This was discussed in detail in I2ll (we also refer 
the reader to the classical articles [13, l43j |) , and the phe- 
nomenology of this melting process is strikingly similar 
to the one of the dynamics of super-cooled liquids. 

In particular the melting time diverges super- 
exponentially as the first order phase transition is ap- 
proached, as can be understood by standard nucleation 
arguments. In glasses the equilibration time also grows 
super-exponentially, and according to some diverges at 
an ideal glass transition temperature. In a system with a 
first order transition we also observed diverging dynam- 
ical and static correlation lengths. In particular, the dy- 
namical length scale is associated to heterogeneities 
in the dynamics. It uses a four-point density correla- 
tor in both time and space, and it led to the notion of 
the so-called dynamical susceptibility (usually refereed 
to as Xi)- The X4 susceptibility was observed to grow 
also in glass formers [1, |9| . There is also a static (ther- 
modynamic) growing length scale associated to a first 
order phase transition, the so-called point-to-set corre- 
lation which is the correlation of a sub-system with its 
frozen boundaries. This correlation length was observed 
to grow also in glass formers (lol - [l^ . 



Despite the profund analogy [21[ there are crucial dif- 
ferences between glassy dynamics and the melting one 
through a standard first order transition. Consider ther- 
modynamic properties: there is no latent heat associated 
to the glass transition, as the energy at a glass transition 
is continuous, while in general there is a latent heat in 
first order phase transitions. This difference can be over- 
come by considering a first order transition driven purely 
by entropy where there is no latent heat, as for instance 
in hard spheres [ill, liquids crystal and other systems 
[4^. A more serious difference between glassy dynamics 
and standard melting appears in the dynamical behav- 
ior. Usually the melting process happens once for all and 
when the system has melted, it stays in the liquid phase, 
usually a very different one from the ordered phase. The 
equilibrium glassy dynamics, on the other hand, is a sta- 
tionary process which is time-translationally invariant. 

Considering the results of the previous sections, one 
sees that both these differences vanish on the Nishimori 
line: the melting dynamics is equivalent to the stationary 
equilibrium dynamics. Moreover, there is no latent heat 
as the energy is given by ([T5| . More importantly, if there 
is a first order phase transition on the Nishimori line 
then there are divergent length and time scales, and the 
mapping between melting and equilibrium dynamics tells 
us that there are genuine divergences in the equilibrium 
dynamics as well! To conclude: if there is a first order 
phase transition on the Nishimori line then we found an 
exact and simple-to-study realization of the ideal glass 
transition. This is the main thesis of this paper and we 
shall now pursue it, first in a mean-field system, and then 
in a finite dimensional one. 
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II. GLASSY DYNAMICS AS A MELTING 
PROCESS IN MEAN FIELD SYSTEMS 

We shall first study the p-spin model on the Nishimori 
line in the mean-field setting. As we shall see, this will be 
equivalent to the usual spin-glass mean field theory, so it 
is useful to first review the properties of mean-field spin 
glasses. We will consider a random lattice where every 
spin is involved in exactly c ~ 0(1) interactions. Such a 
diluted p-spin model is called the XOR-SAT problem in 
the literature and can be solved using the replica or the 
cavity method (see [H, |4^, |47|). We chose to work with 
the diluted p-spin model instead of the more common 
fully connected one because it can be simulated in a 
time linear (versus quadratic) in the size of the system, 
and because a distance between spins is naturally defined. 

The Ising p-spin glass (22l - [2^ , that is Hamiltonian ([1]) 
when p = 1/2, provides a mean-field theory for the struc- 
tural glass transition [l^ . The thermodynamic behavior 
of the p-spin glass undergoes the following changes as the 
temperature T is decreased (see Fig.[T]): At high temper- 
ature, the system is in a paramagnetic (or liquid) phase. 
Below the so-called dynamical glass temperature, that we 
denote here Tmct, the paramagnetic state shatters into 
exponentially many Gibbs states: the energy landscape 
is therefore divided into exponentially many dynamically 
attractive regions, all well separated by extensive free- 
energetic barriers. This leads to a breaking of ergodicity 
on any non-exponential time-scale and to the power-law 
divergence of the equilibration time at Tmct 
This is the equivalent of the mode-coupling transition in 
structural glasses. Note, however, that this dynamical 
transition is not a phase transition in the usual sense: 
there is no non-analyticity in the free energy at the tran- 
sition. It is only a topological transition in the config- 
urational space that affects the dynamics of the system 
(thus the name dynamical transition). 

As the temperature is further lowered, the number 
of states (relevant for the Boltzmann measure) decays. 
A second static Kauzmann transition is then reached at 
Tk when the number of relevant states becomes sub- 
exponential (and in fact finite) and the structural entropy 
(or complexity) vanishes. The RFOT departs from this 
mean-field picture and argues that in finite dimensional 
systems Tmct is only a crossover so that the relaxation 
time diverges at Tk- But in this section we will restrict 
ourselves to the mean-field case. 

Another property of the p-spin model with p — 1/2, 
that will turn out to be very useful in what follows, is 
that the annealed averages are equal to the quenched ones 
above the Kauzmann temperature. 



modynamic average over configurations is taken. Only 
after that the disorder average is taken, and this is the 
physically correct way to proceed in most disordered sys- 
tems. In the annealed case (which is in general only 
approximate) the average over disorder is taken always 
at the same time as the average over configurations. 

In the p-spin model above the Kauzmann transition, 
Tk, the quenched free energy is asymptotically equal to 
the annealed one -/3/ = [log {Z)/N]^/^ = log [Z]^^r^/N + 
o(l), see [H, [l^]. Assuming at least exponentially rare 
large deviations for thermodynamic quantities, we can 
see that choosing first a disorder realisation and then a 
random configuration of a given energy is the same as 
taking a random configuration and choosing the disor- 
der at random such that the configuration has the same 
energy. This is because in the second case instances of 
disorder are chosen proportionally to the value of the par- 
tition function at a corresponding temperature but since 
[log (^)]i/2 — -^OS [^]i/2 + o(iV) this second way chooses 
typical instances as wellQ. 

In what follows in the mean field p-spin model we can 
thus freely exchange the quenched average of thermody- 
namic quantities for annealed one or vice versa, i.e. 



(27) 



We stress that the above holds only for disorder aver- 
ages of thermodynamic quantities (i.e. when the quan- 
tity or its density is bounded by an A^-independent con- 
stant), for instance it does not have to hold that [Z]^^^ — 
[Z^] ^^2- We now denote the number of interaction by M 
and compute explicitly the free energy. First, we obtain 



2M 



2^ cosh*'^ (/3) . (28) 



The free energy density for T > Tk is thus simply 

M 

-/3/ = log2 + — log(cosh/3), 

and the energy is given by 

d(/3/) _ M 



d/3 N 



— tanh(/3). 



(29) 



(30) 



as we could also obtain from eq. ((27)) . These results are 
correct for all T > Tk- Just to be concrete, for the 3-spin 
model where every spin is involved in c = 5 interactions, 
one finds Tk — 0.849, while the system is in the many 
valleys glassy phase for all Tk <T < Tmct = 0.936. 



A. Annealed and quenched averages in the mean 
field p-spin model 



B. The Nishimori line in a mean field p-spin model 



In the quenched average the disorder realization (i.e. 
the lattice and signs of interactions) is fixed, and the ther- 



Let us now consider the general identity ([9]) , when the 
quenched and annealed averages at inverse temperature 
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Kp are equal we arrive to 

[A{{J^,k)AS^},P)]p = [{A{{J^JkCT^a,ak}dS^'y^},P))K^. 

(31)' 

where {■)k is the thermal average at inverse temper- 
ature Kp. This means that physics of the mean- field 
spin glass at temperatures /3 — Kp is uniquely mapped 
on physics of the same model but with a ferromagnetic 
bias pnl corresponding to the Nishimori line. This holds 
above the Kauzmann temperature, T > Tk, and in any 
model where the annealed average equals the quenched 
one (that does not include any finite-dimensional model 
we are aware of). 

We give several more specific examples of the above 
statement. According to the previous section the melting 
process on the Nishimori line is equivalent to the equilib- 
rium dynamics on the Nishimori line. In the mean field 
p-spin moreover it is equivalent also to the equilibrium 
dynamics in the spin glass problem, p = 1/2, as long as 
T>Tk. One has 



Pnl(;3') 



.it) = Cp^dt) 



f 



J-PNL(^i') 

Jfp 



(/3,/3',TO) 



(32) 
(33) 



) = /^i?(/3,/3',™). 



(34) 



The above mapping allows us to understand in a simple 
manner the glassy behavior of the p-spin model 75| . Let 
us consider an equilibrium configuration at temperature 
T > Tk, the time needed to decorrelate diverges at the 
dynamical, or mode-coupling transition Tmct- The cor- 
relation function develops a plateau whose length 
diverges at Tmct- Using the above mapping, these fea- 
tures translate in the system on the Nishimori line. On 
the Nishimori line there is a a first order phase transi- 
tion at Tp — Tk from a paramagnetic to a ferromagnetic 
phase with the following properties. Initializing the dy- 
namics in the fully ordered configuration, the magnetiza- 
tion decays to zero only for temperatures larger than the 
spinodal one Tg — Tmct (since for T < Tg the system 
will be trapped in the ferromagnetic state). The time 
needed to relax to zero magnetization diverges at Tg as 
a power law (since we deal with a mean-field system). 
The magnetization m(t) develops a plateau whose length 
diverges at Tg. 



C. Simulating glassy dynamics as a melting process 

We shall now simulate the melting process on the 
Nishimori line. A useful application of the idea discussed 
above is that it provides efficient numerical simulations 
of the equilibrium dynamics in the spin glass. Instead 
of equilibrating Monte-Carlo simulation in the spin glass 
model we can simply initialize in the fully magnetized 
configuration on the Nishimori line. 

We have performed Monte-Carlo simulations of the 
melting process of a large system on the Nishimori line. 



For the 3-spin model with c = 5 for temperature /?, we 
have prepared a random realization of the system with 
N = 90000 spins and a proportion of negative couplings 
given by the Nishimori condition eq. ([3]), and initialized 
the simulation with all spins Si — 1. We then use Monte- 
Carlo dynamics with the Metropolis update rule. The 
results are shown in Fig. [2] The exact mapping ensures 
that the magnetization evolves the same way as the cor- 
relation function in the spin glass, and that the dynamic 
susceptibility evolves as the 4-point correlation function 
in the spin glass. The relaxation time diverges as 



r cx 1 



T 



MOT 



-2.8(2) 



(35) 



where the critical exponent was obtained by fitting the 
curve in the inset of Fig. [2] and agrees well with the val- 
ues obtained in [Tl| . The maximum of the dynamical 
susceptibility diverges as 



1 - 



T 



MCT 



-1.15(15) 



(36) 



where the critical exponent was obtained by fitting the 
curve in the inset of Fig. [21 This agrees quite well with 
the predictions of the mode-coupling theory d, ISTt - fsoj 
where the exponent is equal to —1. 

Let us stress that this simulation would take a much 
longer time without the above trick that gives us an equi- 
librium configuration for free. Moreover, the method ex- 
tends to temperatures Tk < T < Tmct where Monte- 
Carlo equilibration takes an exponential time and is 
hence infeasible for such large system sizes [76]. 

The above mapping also allows to understand the 
point-to-set correlations [ll|) in a simpler way. Point- 
to-set correlations are defined as follows: freeze the sys- 
tem in an equilibrium configuration. Then consider a 
large compact droplet in the system and un-freeze it to 
observe the infiuence of the frozen boundaries. For each 
temperature, one defines the equilibrium point-to-set vol- 
ume as the largest size of the droplet such that the bound- 
aries are correlated with the center of the droplet. Us- 
ing the above mapping, this notion translates to a much 
simpler, and more easy to simulate, one: we just have 
to consider a system with fully ferromagnetic boundary 
conditions and ask if the magnetization stays positive in 
a center of a droplet of a certain size: This is the very 
same question we addressed in the purely ferromagnetic 
system in [2ll |. 

Point-to-set correlations in the mean-field spin glass 
case can be derived analytically and obey a simple solv- 
able recursion (see [ll|). Our mapping shows (as can also 
be seen explicitly in the equations) that this recursion is 
identical to the one describing the way the magnetiza- 
tion decays at distance d from ferromagnetic boundary 
conditions on the Nishimori line. We plot the solution to 
these equations in Fig. [31 Actually since the equations in 
question are very similar to what happens in pure ferro- 
magnets, it is not surprising that the decorrelation length 
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FIG. 3: Magnetization at distance d from a ferromagnetic 
boundary conditions in 3-spin model on a Bethe lattice of 
coordination c = 5 on the Nishimori line. This is the point- 
to-set correlation in the spin glass case. The correlation length 
diverges with exponent 1/2 at Tmct = 0.936. 



FIG. 2: Dynamic ferromagnetic susceptibility and decay of 
the magnetization m(t) in the melting on the Nishimori line 
for the 3-spin model on the Bethe lattice with c = 5 and TV = 
90000 spins, averaged over 500 realizations of the dynamics 
except for the three lower temperature where we used only 50 
runs. From left to right; T = 1.3, 1.2, 1.1, 1.05, 1.0, 0.99, 0.98, 
0.97, 0.96, 0.955 and 0.95. Bottom inset: the melting time 
[defined by m(r) = 1/e] that diverges as a power law p5|l . 
Top inset: The maximum of the ferromagnetic susceptibility 
X F diverges as (|36|l . The plotted quantities are identical to the 
dynamical susceptibility X4(i), eq. p3)) . and the equilibrium 
autocorrelation C{t), eq. (I32|l . in the spin glass phase. 



grows when Tmct is approached with the exp onent 1/2 
that was first obtained, for this model, in [ll|. 



III. GLASSY DYNAMICS AS A MELTING 
PROCESS IN A FINITE DIMENSIONAL MODEL 

We now discuss how is this mean field picture modified 
in finite dimensional systems. Since we are dealing with 
a glass transition, we should in principle use the mosaic 
theory, which is the counter-part of nucleation theory 
for glasses 

[i[i3,[il-[I3- The fundamental point in this 
work is that due to the mapping to melting above a first 
order phase transition (from a ferromagnetic phase to a 
paramagnetic one) it is sufficient to consider the standard 
nucleation theory. This is a great conceptual simplifica- 
tion with respect to the many studied lattice models of 
glass formers [33l - l35| . 



A. Nucleation theory for the melting of a 
disordered magnet 

Let us very briefiy recall the basics of nucleation theory, 
and discuss the effect of disorder for an entropy-driven 
first order phase transition. Nucleation stems from a 
competition between the bulk entropy difference and the 
surface tension between two phases. When a system is 
in a metastable phase the total free-energy cost /S.F{1) 
of a compact droplet of size i is given by a combina- 
tion of the bulk entropy —i'^TSs and the surface tension 
£^T, where Ss is the entropic difference between the two 
phases, < d — 1 is an exponent characterizing the free 
energy cost of the surface, and T is the surface tension 
(we omit constant pre- factors for convenience) . In a non- 
disordered systems we have 9 = d—1, but when disorder 
is present the shape of the droplet can adapt in order to 
take advantage of the local configuration of the disorder, 
so for the best droplet of each size we have 9 < d—1. The 
two contributions balance when the droplet size is 

/ r \ ^ 

4 cx i^—j cx (T - TJ-'^ (37) 

with 1/ — l/{d~9), where we assumed that the differ- 
ence ds is linear when approaching the first order phase 
transition. For sizes larger than £c fiipping a droplet is 
advantageous for the system. In order to flip, the ther- 
mal fluctuations must cross the free energetic barrier. In 
a non-disordered system one would expect that the bar- 
riers grow as i^. However, this is only the free energy 
of the best droplet along the way, but the path that the 
dynamics is following is not bound to pass through this 
best droplet; it might actually have to go through larger 
barriers, and it is thus safer to assume instead that the 
barriers grow as with ip > 9. Using the Arrhenius 
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formula we thus obtain the nucleation time 



T oc exp 



A 



exp 



A 



(T - ^,)^ 



(38) 



with 7 = ijj/{d - 0). 

So far we have, however, neglected an important in- 
gredient: disorder induced fluctuations. We shall now 
present a scaling analysis in order to estimates the ef- 
fect of the disorder on a first-order transition. Consider 
a droplet of size £, we expect that it can gain a free en- 
ergy of order AFdis — SJ£'^^'^ just from fluctuations in 
the free energy due to the disorder variance SJ. The 
point is that this might be enough to make the droplet 
flip if the disorder induced term is larger than the sur- 
face term Ti^ . Indeed, if 6* < d/2 there is a finite length 
^dis cx (r/(5J)^^^'^^^ beyond which the droplet flips in- 
dependently of temperature. Hence in such a case the 
critical droplet size (that is, the point-to-set size) does 
not diverge at the transition point, and consequently the 
transition cannot be of first order. If we want a genuine 
first order phase transition, we thus must have 9 > d/2, 
or equivalently 



ly = l/{d~0) > 2/d. 



(39) 



For instance, the introduction of small amount of disorder 
in the pure ferromagnetic model, where 9 — d~l, cannot 
lead to a first order phase transition for d < 2. The sit- 
uation is the same in the marginal case d = 2, as proven 
rigorously by Aizenmann and Wehr [39l |. According to 
this analysis, we thus need at least a three-dimensional 
model if we want to have a hope for a genuine first order 
phase transition in the presence of disorder. 

Finally, following the discussion from a companion pa- 
per [131 , let us emphasize that the relevant mechanism for 
melting is nucleation and growth. The nucleation time 
T ([55]) only gives the average time needed to nucleate a 
critical droplet at a given position in the system. If the 
system is much larger than the critical nucleus size the 
melting process is not done via the growth of a single 
droplet but instead results from the nucleation of many 
droplets that will simultaneously grow and invade the 
system. This growth time has to be added to the nucle- 
ation time to obtain the relaxation time needed to exit 
from the metastable state in a very large system. The 
probability to nucleate a critical droplet by a unit of vol- 
ume and time is proportional to e~^^^(''=) (x 1/t. Hence 
in a very large system there is roughly one droplet per 
volume T after a unit of time. Without disorder, the vol- 
ume of large droplets is expected to grow algebraically 
with time and consequently in the systems without dis- 
order the total relaxation time is of the same order as the 
nucleation time r, and the typical size of grown droplets 
when they percolate is also scaling as t. This can be 
observed by looking to the dynamical ferromagnetic sus- 
ceptibility, defined in eq. as we have discussed in 
detail in [21] for the Potts model. Clearly, the moment 
when such droplets percolate denotes the moment when 



the system is most heterogeneous, where roughly half of 
it is in the liquid phase and the other part still in the or- 
dered one, so that the two-point dynamical ferromagnetic 
susceptibility is maximal. 

Disorder will modify the growth behavior. Instead of 
the algebraic nucleus growth we expect instead an acti- 
vated growth due to pining of the interfaces, just like in 
the coarsening of disordered magnets where the growth is 
logarithmic in time [g^l • The interplay between activated 
growth and the appearance of new droplets makes the to- 
tal analysis more involved than in the case without disor- 
der, and is a subject worthy of closer inspection in order 
to understand the melting of a crystalline ordered phase 
in presence of disorder. If the growth is slow, the melt- 
ing will have to wait until many droplets have nucleated 
everywhere. Again, this gives a time of order r for the 
melting process, but when droplets percolate, they are 
typically much smaller than in the non-disordered case. 
We will see that the dynamical magnetic susceptibility is 
again a suitable tool to quantify this understanding. 



B. Nucleation vs mosaic on the Nishimori line 

We now focus on a first order phase transition on the 
Nishimori line. Since we are discussing an entropy driven 
transition in a system with disorder, the analysis of the 
former paragraph applies. On the Nishimori line, how- 
ever, the melting process is equivalent to the equilibrium 
dynamics. That means the system that is in a given 
equilibrium configuration decorrelates from this config- 
uration by the nucleation mechanism to reach another 
equilibrium configuration which again will be left by ac- 
tivated nucleation and so on. Such a stationary melting 
process, where one transits from one equilibrium con- 
figuration to another by nucleation is nothing else than 
the phenomenology of the RFOT and the mosaic pic- 
ture. Due to our mapping, however, it simply appears as 
a consequence of standard nucleation arguments. 

For instance, e g . ( [ 551) is nothing but the Adam-Gibbs 
relation [l^, [isl - llTj . while eq. ([57]) is the mosaic length. 
The entropy difference is now associated with the number 
of metastable state so that Ss = T, (where E is the config- 
urational entropy). According to the condition p9p . we 
need > 2/d, a conclusion was also reached in the sem- 
inal article on the RFOT 17] where in fact it was even 
argued that v = 2/d, which combined with ^ = 9 gives 
7 = 1, which gives the original Adam-Gibbs relation. If 
indeed i/ = 2/d it suggests that if there is a first-order 
transition on the Nishimori line, it is a marginally stable 
one with respect to the smoothing by disorder. These 
questions, and in particular whether indeed v = 2/d ot 
not, deserve further investigations. 

An interesting (and hopefully clarifying) relation in our 
mapping is given by eq. (|27p . The four-point dynami- 
cal susceptibility Xi{t)j that describes the heterogeneous 
glassy dynamics, is simply the dynamical ferromagnetic 
susceptibility in the melting process. Xii^) thus acquires 



10 



a simpler interpretation, as discussed in section IIII Al 
The mapping on the Nishimori hne thus allows to dis- 
cuss the dynamics of glass forming liquids as a much 
simpler melting process. Let us illustrate this in a three- 
dimensional model. 



C. A finite dimensional p-spin model 

There are many ways how to introduce a p-spin model 
on a finite dimensional lattice (see for instance [6l|). We 
will use the following model on a three-dimensional cubic 
grid with each spin being involved in 10 different 5-body 
interactions 

^ ^ t/j iSi^up^lcft bright "Sback 

i 

"^i 'S'iS'down'S'loft'S'i-ight'S'front J (40) 

i 

where the indexes indicate the position of the spin on 
the grid with respect to the spin i. In the mean field case 
— that is on a Bethe lattice with p — 5 and c — 10 — this 
model can be solved using the cavity method and it has 
qualitatively the phase diagram from Fig. [1] with a first 
order phase transition at Tc (p = 1) = 2.4 in pure case 



while on the Nishimori line, 
Tc(pnl) = 0.96. 



=2/3 



p/{l — p), we have 



D. Numerical results 

We now present results of numerical simulations of 
Hamiltonian (j40p on the three-dimensional cubic grid on 
the Nishimori line. Studying the equilibrium dynamics 
on the Nishimori line, using the results of section U is 
quite simple. We simulate the melting process starting 
from the fully ordered configuration, just as we did in 
the mean field system. Of course now the system on 
the Nishimori line is not equivalent anymore to the spin 
glass case with p = 1/2, but this is irrelevant for our dis- 
cussion. If there is a first order phase transition on the 
Nishimori line, then the melting process with a Vogel- 
Fulcher-Talman-like relaxation time and diverging length 
scales implies that the equilibrium dynamics undergoes 
an ideal glass transition described by the RFOT. 

The correlation function for equilibrium dynamics (or 
equivalently the magnetization in the melting process) is 
plotted versus temperature in Fig. |4l It yields the stan- 
dard picture with a growing plateau and the melting time 
— or equivalently the equilibrium auto-correlation time — 
grows faster than exponentially, and is perfectly consis- 
tent with a VFT divergence at Tc — 0.88, as attested in 
Fig. [5] Note that we use the usual form of the VFT, but 
we could have used the generalized form eq. ([55]) with 
similar agreement but with slightly different values of Tc- 

The upper part of Fig. |4] shows the dynamical suscep- 
tibility or equivalently the 4-point dynamical susceptibil- 




FIG. 4: The dynamics on the Nishimori line for the 3- 
dimensional 5-spin model. Top: Dynamical susceptibility or 
equivalently the 4-point correlation function, X4(^) — Xf (^)i 
as a function of time. Bottom: The magnetization or equiva- 
lently the equilibrium autocorrelation function, C{t) = m{t), 
as a function of time. We use 1000 realizations of the dynam- 
ics, except for the two lowest temperatures where we used 
only 50 realizations. From left to right: T = 2.2, 2.0, 1.8, 1.6, 
1.5, 1.4, 1.3, 1.25, 1.2, 1.15, 1.1, 1.075, 1.05 and T = 1.025. 
The peak of the dynamical susceptibility seems to increase 
slower than a power law of the relaxation time r. We show a 
power-law line with exponent 0.22, which is the best power- 
law estimate for the largest time; however such a low value, 
and the general trend of the maxima, suggest that asymptot- 
ically the growth is logarithmic in r (and thus algebraic in 
T — Tc) as expected for an actived coarsenning. 



ity. Its maximum is growing with the growing re- 



laxation time r. In [21| we observed for the ferromagnetic 
two-dimensional non-disordered case that xT'^^ ~ The 
present data for the 5-spin model on the Nishimori line 
instead indicate a low-exponent power law; in fact for 
the largest time the best power-law fit gives xT'^^ w r°-^. 
Such a low exponent is compatible with a logarithmic 
growth in time; in fact a similar behavior is observed in 
the activated coarsening of the random field Ising model 
and in the random bond Ising model, where a logarith- 
mic growth is expected 60] but in numerical simulations 
a power-law with low exponent is observed [gl] . The low- 
exponent power law is in this case usually attributed to 
pre- asymptotic effects. The slower divergence of xT^^ 
the presence of disorder (to be compared with the non- 
disordered transition in the Potts model observed in |21|) 
is thus due to an activated coarsening for the growth of 
the nuclei, as expected. 

We also investigated the point-to-set correlation func- 
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FIG. 5: Melting time or equivalently the relaxation time r, 
in the 3-dimensional 5-spin model as a function of the tem- 
perature T (plotted in a log-linear way on the left, and loga- 
rithm of the time versus 1/{T ~Tc) on the right) . The depen- 
dence is very well fitted with a Vogel-Fulcher-Talman form 

T oc e(^-^<=) , with Tc « 0.88. 
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FIG. 7: Point-to-set correlation length, defined as the size of 
the box where the magnetization in Fig. [6] goes below 2/3. 
It is well fitted by 1/{T - Tc)" with ly 0.75 and Tc = 0.9. 
Using instead the 1/[T{T - Tc)]" from (O gives ly ^ 0.67. 




the 5-spin model on the S-dimensional grid on the Nishi- 
mori line support the existence of a first order transition 
close to Tc w 0.9. Since the melting process is equivalent 
to equilibrium dynamics on the Nishimori line, this is a 
perfect realization of RFOT ideas in a finite dimensional 
setting. Below the Kauzmann transition the equilibrium 
physics is, however, given by the ferromagnetic phase, 
not by an ideal glass phase. This is not a problem, as 
this is after all what happens in real systems — the crys- 
tal is always the correct equilibrium state — and as we 
are mainly interested in how a liquid becomes a glass, 
not in the out-of-reach glass state below Tk (which is 
anyway ill- defined j35|). 



FIG. 6: Point-to-set correlation, or equivalently equilibrium 
magnetization m{T, L) in a box of size with ordered bound- 
ary conditions, for different temperatures. A clear sign a 
growing length scale as the temperature is lowered is observed. 



tion [13|. We consider a box of size N = L^, apply 
fully ferromagnetic boundary conditions, and measure 
the equilibrium magnetization inside the box for different 
temperatures and sizes. Again, since we are on the Nishi- 
mori line, this is nothing but the point-to-set equilibrium 
correlation length, computed in a much faster way. The 
data are plotted in Fig. \6\ and the correlation length in 
Fig. [T] We observe a growing length scale compatible 
with a power law divergence at Tc- By fitting our data 
using the formula e oc [T{T - Tc)]'" of eq. ^ we found 
that Tc is compatible with the value 0.88, and that the 
exponent v is compatible with the value 2/3 advocated 
in [13]. 

Given we barely span two decades in Fig |7] one cannot 
be truly conclusive. Still, all our numerical results for 



E. Is there really a first order transition? 

The numerical data are compatible with the existence 
of a first order phase transition on the Nishimori line, 
and as a consequence, with the existence of the ideal 
glass transition in a finite dimensional system. But is it 
really a genuine first order phase transition? Of course, 
this cannot be answered based only on numerical simu- 
lations, and we cannot exclude that if one goes to much 
larger times and to point-to-set lengths equal to, say, 
^ « 100 a crossover to a second-order transition arises. 
Or worse, there might be no phase transition at all. Note, 
however, that such first-order transitions definitely exist 
at the mean field level (see section Hi B[) and that at least 
second-order ones do exist on the Nishimori line (see [63|) 
in finite dimension: We are thus assured that a transi- 
tion can be found. At least, the issues and criticisms 
discussed in [g^l, where it was argued that the RFOT 
transition in finite dimension maps only to a spin glass 
droplet model in a field, where it seems there is no tran- 
sition at all ^Bj, do not apply in our case. The difficulty 
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here is, however, to decide if on the Nishimori hne there 
can be a first order phase transition in finite dimension. 

Unfortunately, in all the solvable cases we found in 
the literature the transition on the Nishimori line is of 
the second order, e.g. [6^. This is not very surprising 
because all these cases concern either 2-dimensional sys- 
tems to which the theorem of Aizenman-Wehr about 
smoothening of first order phase transition in the pres- 
ence of disorder applies, either system which have a 2nd 
order phase transition even in the absence of disorder. 
If a genuine first order transition can be found on the 
Nishimori line in a 3D system is hence a challenging open 
question, which is clearly worth investigating. 

Proving the existence or impossibility of a first order 
phase transition in a finite dimensional model on the 
Nishimori line would be a nice achievement. A positive 
answer would prove that a true VFT divergence exists in 
finite dimension, and that it is associated with a RFOT 
mechanism. However, on the practical side, this is not 
such a fundamental issue. If the transition becomes of 
2nd order at sizes so large that the equilibrium behavior 
cannot be observed in a time smaller than any experi- 
mental time, then the first order approach for the melting 
problem (and the RFOT one for the glassy dynamics) is 
the correct description at the observable times scales. 



IV. DISCUSSION 

In this paper we have further studied the analogies 
between the melting dynamics above a first order phase 
transition and the equilibrium glassy dynamics discussed 
in [2l| by exploring a class of models where the two pro- 
cesses are strictly equivalent: Ising spin systems on the 
Nishimori line (some disordered Potts models share the 
properties discussed in this paper \63\ ) . In the light of 
the mapping between melting and equilibrium dynamics, 
we investigated several features of the dynamics of super- 
cooled liquids: the diverging relaxation time, the plateau 
in the correlation function, dynamical heterogeneities re- 
lated to the dynamical susceptibility X4, and the point- 
to-set correlations. We show that they all can be easily 
recovered, understood, and studied in the much simpler, 
and more familiar, setting of first order phase transition. 
Since the RFOT was invented in analogy with first or- 
der transitions, it is rather interesting that there exist 
systems where the two are equivalent. 

In particular, on the Nishimori line, the dynamical — 
or mode-coupling — transition corresponds to the spin- 
odal point, while the Kauzmann transition corresponds 
to the first order ferromagnetic transition. The mosaic 
approach is replaced by simpler nucleation arguments. 
We hope this approach will help in expanding our com- 
prehension of the glass transition problem. To conclude 
this work, we would like to make several comments. 

Ideal glass transition? A very fundamental open 
question in the field of glassy systems is the existence of a 
finite dimensional model where one can observe the VFT 



divergence. The bottom-line of our result here is that if 
there is a first order phase transition on the Nishimori line 
then the ideal glass transition exists. It would be very 
interesting to prove (or disprove) the existence of such 
a transition on the Nishimori line. Hopefully, this will 
be easier than the original question, as on the Nishimori 
line we are discussing a simpler to study first order phase 
transition. Moreover, the gauge transformation can be 
used to prove many results rigorously j26j . 

String-like motion It is well known in the theory 
of nucleation and first order phase transitions that the 
nucleating droplets can also be non-compact and fractal 
[4^ . As noted by W], the stringy nuclei observed close 
to the pseudo-spinodal \§7} in first-order transitions could 
be the analog of the string-like motion observed in the 
relaxation in glassy systems The results presented 
here and the exact analogy between melting and glasses 
thus provide a firmer motivation for this. It would also 
be interesting to observe the string-like motion during 
the melting process of superheated solids: This is maybe 
related to the rings and loops observed in bulk melting 
simulations [69| . 

Mode coupling transition as a mean field ap- 
proach Many theoretical results about glasses are 
based on the mode coupHng theory [13, [11]. The MCT 
predictions are exact in the mean-field p-spin model, 
and in this model the mode coupling transition can be 
mapped to a spinodal point of a first order phase tran- 
sition on the Nishimori line. One hence cannot escape 
the conclusion that the transition described by MCT is a 
kind of mean-field melting, and the mode coupling tem- 
perature is a corresponding spinodal point. Of course, 
this suggestion is present in the literature since the work 
of 14-16], but our mapping makes it very explicit. 

Simulating MCT behavior A straightforward but 
powerful application of our mapping is to perform fast 
simulations of the equilibrium glassy dynamics in mean- 
field models, as we have done for the p-spin, by simply 
starting from the fully ordered configuration. This should 
allow us to go beyond the current size limitations and 
investigate better numerically the mean-field systems and 
the finite-size corrections (in the spirit of e.g. ^7^). This 
is in fact what we have started to do in section Hi CI 

Crossover MCT/finite dimension One of the 

crucial issues in the RFOT is to understand better nucle- 
ations in glassy systems, and our mapping allows to put 
this question back into the standard framework of first 
order phase transitions. The crossover from the mode 
coupling behavior to the finite dimensional activated dy- 
namics (described e.g. by the mosaic theory) then be- 
comes the usual crossover from spinodal to activation at 
a first order phase transition ^42.] . A promising direction 
in which these ideas could be extended is the Kac limit 
[tH , or the field theoretic computation of instanton in 
order to estimate the free energetic barriers. Note also 
that the system is guarantied to be replica symmetric on 
the Nishimori line which may simplify many calculations. 



13 



Acknowledgments 

It is a pleasure to thank G. Biroli, J-P. Bouchaud, A. 
Cavagna, S. Franz, J. Kurchan, J. Langer and H. Yoshino 



[1] P. G. Debenedetti and F. H Stillinger, Nature 410, 6825 

pp. 259-267 (2001). 
[2] M. D. Ediger, C. A. Angell, S. R. Nagel, J. Phys. Chem. 

100, 13200 (1996). 
[3] C. A. Angell, Journal of Physics and Chemistry of Solids 

49, n. 8, 863-871 (1988). 
[4] H. Vogel, Phys. Z. 22, 645 (1921). G. S. Fulcher, J. Am. 

Ceram. Soc. 8, 339 (1925). 
[5] W. Kauzmann, Chem. Rev., 43 (2), 219 (1948). 
[6] G. Adam, J. H. Gibbs, J. Chem. Phys. 43, 139146 (1965). 
[7] S. Franz and G. Parisi, J. Phys.: Condens. Matter 12, 

6335 (2000). 

[8] C. TonineUi et al, Phys. Rev. E 71, 041505 (2005). 
[9] E. R. Weeks et al, Science 287, 627-631 (2000). 
[10] J.-P. Bouchaud and G. BiroU, J. Chem. Phys. 121, 

73477354 (2004). 
[11] A. Montanari and G. Semerjian, J. Stat. Phys. 124, 103 
(2006). 

[12] A. Montanari and G. Semerjian, J. Stat. Phys. 125, 23 
(2006). 

[13] A. Cavagna, T. S. Grigera and P. Verrocchio, Phys. Rev. 
Lett. 98, 187801 (2007). G. Biroh, J.-P. Bouchaud, A. 
Cavagna, T. S. Grigera and P. Verrocchio, Nature Physics 
4, 771 (2008). 

[14] T. R. Kirkpatrick and D. Thirumalai, Phys. Rev. Lett. 

58, 2091 (1987). T. R. Kirkpatrick and D. Thirumalai, 

Phys. Rev. B 36, 5388 (1987). 
[15] T. R. Kirkpatrick and P. G. Wolynes, Phys. Rev. A 35 

3072 (1987). 

[16] T. R. Kirkpatrick and P. G. Wolynes, Phys. Rev. B 36 
8552 (1987). 

[17] T. R. Kirkpatrick, D. Thirumalai, P. G. Wolynes, Phys. 

Rev. A 40, 10451054 (1989). 
[18] M. Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999). 

[19] J.S. Langer, Phys. Today 89 (2007)^ 

[20] G. Biroh and J.-P Bouchaud, arXiv :0912,2"542] (2009). 
[21] F. Krzakala and L. Zdeborova, Glassy aspetcs of melting 

dynamics, in preparation. 
[22] B. Derrida, Phys. Rev. Lett. 45, 79 - 82 (1980). B. Der- 

rida, Phys. Rev. B 24, 2613 - 2626 (1981). 
[23] D. J. Gross and M. Mezard, Nucl. Phys. B 240, 431 

(1984). 

[24] F. Ricci-Tersenghi, M. Weigt and R. Zecchina, Phys. Rev. 
E 63, 026702 (2001). S. Franz et. al. Europhys. Lett. 55, 
465-471 (2001). 

[25] H. Nishimori, Prog. Theor. Phys. 66, 1169 (1981). 

[26] H. Nishimori, Statistical Physics of Spin Glasses and In- 
formation Processing: An Introduction (Oxford, 2001). 

[27] Y. Ozeki, J. Phys. A: Math. Gen. 28, 3645 (1995). 

[28] Y. Ozeki, J. Phys.: Condens. Matter 9, 11171 (1997). 

[29] S. F. Edwards and P. W. Anderson 1975 J. Phys. F: Met. 
Phys. 5 965 

[30] M. Mezard, G. Parisi, and M. A. Virasoro, Spin-Glass 
Theory and Beyond, Lecture Notes in Physics Vol. 9 
(World Scientific, Singapore, 1987). 



for interesting discussions and comments. 



[31] J. Hertz and K. H. Fischer, Spin Glasses (Cambridge 

Universtity Press, 1991). 
[32] M. Goldstein, J. Chem. Phys. 51, 3728 (1969). 
[33] A. Lipowski, and D. Johnston, Phys. Rev. E 64, 041605 

(2001) 

[34] M. R. Swift, H. Bokil, R. D. M. Travasso and A. J. Bray, 

Phys. Rev. B 62, 1149411498 (2000). 
[35] A. Cavagna, I. Giardina and T. S. Grigera J. Chem. Phys. 

118, 6974 (2003). 
[36] A. Georges, D. Hansel, P. Le Doussal and J-M Maillard, 

J. Phys. 48, 1 (1987). Y. Kasai and A Okiji, Prog. Theor. 

Phys. 69, 20 (1983). 
[37] S. Franz and M. Leone, J. Stat. Phys. Ill, 535 (20903). 
[38] F. Guerra and F. L. TonineUi, Commun. Math. Phys. 

230, 7179 (2002). 
[39] J. Wehr and M. Aizenman, J. Stat. Phys. 60, 287 (1990). 

M. Aizenman and J. Wehr, Phys. Rev. Lett. 62, 25032506 

(1989) 

[40] S. Franz, G. Parisi, Phys. Rev. Lett. 79, 2486 (1997). 
[41] S. Franz et al, arXiv:1001.1746 and arXiv:10 08. 09961 
[42] K. Binder, Rep. Prog. Phys. 50 783-859 (1987). 
[43] K. Binder, Phys. Rev. B 8, 34233438 (1973). 
[44] W. G. Hoover and F. H. Ree, J. Chem. Phys. 49, 3609 
(1968). 

[45] D. Frenkel, Physica A 263 26 (1999)^ 

[46] F. Krzakala and L. Zdeborova, arXiv:0909.3820l to ap- 
pear in EPL (2010). 
[47] L. Zdeborova and F. Krzakala, Phys. Rev. B 81, 224205 
(2010). 

[48] L. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 
(1993). 

[49] J.-P. Bouchaud et al, Physica A 226 243 (1996). J.-P. 
Bouchaud et al, in Spin Glasses and Random Fields, 
edited by A. P. Young (World Scientific, Singapore, 
1998). 

[50] A. Cavagna, I. Giardina and G. Parisi, J. Phys. A: Math. 

Gen. 34 5317-5326 (2001). 

[51] D. Achlioptas and A. Coja-Oghlan, arXiv:0803.2122l 

(2008). 

[52] L Zdeborova and F. Krzakala. larXiv:0902.4185V l (2009). 
[53] F. Krzakala and L Zdeborova, Phys. Rev. Lett. 102, 

238701 (2009). 
[54] S. Franz, G. Parisi, Physica A 261, 317 (1998). 
[55] P. Gillin, H. Nishimori and D. Sherrington, J. Phys. A: 

Math. Gen. 34 2949 (2001). 
[56] H. Nishimori and K. Y. M. Wong, Phys. Rev. E 60, 132 

(1999). 

[57] W Gotze, J. Phys.: Condens. Matter 2 SA201 (1990), 
[58] S. P. Das, Rev. Mod. Phys. 76, 785 (2004). 
[59] G. Biroh and J.-P. Bouchaud, Europhys. Lett. 67, 21 
(2004). 

[60] D. A. Huse and C. L. Henley, Phys. Rev. Lett. 54 (1985), 
2708. D.S. Fisher, D.A. Huse Phys. Rev. B 38 373 (1988). 

[61] D. Alvarez, S. Franz, F. Ritort, Phys. Rev. B 54, 
9756(1996). H. Rieger, Physica A 184, 279 (1992), J. 



14 



Kisker, H. Rieger, H. Schreckenberg, J. Phys. A 27, L853 
(1994). S. Franz and G. Parisi, Eur. Phys. J. B 8, 417422 
(1999). 

[62] R. Paul, S. Puri, and H. Rieger, Europhys. Lett. 68, 881 
(2004). H. Rieger, G. Schehr, R. Paul, Prog. Theor. Phys. 
Suppl. 157, 111 (2005) 

[63] J. Jacobsen and M. Picco, Phys. Rev. E 65, 026113 
(2002). 

[64] M. A. Moore and B. Drossel, Phys. Rev. Lett. 89, 217202 
(2002). 

[65] T. Jorg, H. G. Katzgraber and F. Krzakala, Phys. Rev. 

Lett. 100, 197202 (2008) . 
[66] J. D. Stevenson, J. Schmalian and P. G. Wolynes, Nature 

Physics 2, 268-274 (2006). 
[67] H. Wang, H. Gould and W. Klein, Phys. Rev E 76, 

031604 (2007) 

[68] C. Donati et al, Phys. Rev. Lett. 80, 23382341 (1998). 
[69] X-M Bai and M. Li, Phys. Rev. B 77, 134109 (2008). 
[70] T. Sarlat et al, J. Stat. Mech. P08014 (2009). 
[71] S. Franz and F. L. Toninelli, Phys. Rev. Lett. 92, 030602 

(2004). S. Franz et al. J. Phys. A: Math. Theor. 40 F251- 

F257 (2007). 



[72] This was observed in [3g| in a particular case of eq. @ 
when A = Z, leading to [Z"U^ = [Z"+\^2/[Zh/2- 

[73] As noted by [4l|| there are other possible definitions of 
the dynamical susceptibilities, depending on the order of 
the averages, that lead to slightly different equalities on 
the Nishimori line. For instance, one could use: xpit) ~ 

[{nim] - [{m{t)}f = [{cm] - [{c{t)}f = X4{t). 

[74] Notice that equality of the quenched and annealed av- 
erage holds for a larger class of problems than just the 
p-spin model, see for instance the models where a quiet 
planting is possible [sTI - ls^ . 

[75] It also simplifies many analytical computations. The FP 
potential computed in the spherical p-spin 54] can for 
instance be obtained by considering the free energy with 
a ferromagnetic bias [53]- In the case of Ising spins the 
solution can be obtained readily by looking to the free 
energy on the Nishimori li ne JSol . This was discussed by 
the present authors in [4 q. |47l ]. 

[76] Note that the authors of [ll[ used a similar trick by sim- 
ulating the annealed model. 



